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Abstract 

The instability in the cosmic-ray (CR) precursor of a SN shock is studied. 
The level of turbulence in this region determines the maximum energy of accel- 
erated CRs. The consideration is not limited by the case of weak turbulence. It 
is assumed that Kolmogorov type nonlinear wave interactions together with the 
ion-neutral collisions restrict the amplitude of random magnetic field. As a re- 
sult, the maximum energy of accelerated particles strongly depends on the age of 
a SNR. It can be as high as 10^'^ Z eV in young SNR and falls down to about 10^^ Z 
eV at the end of Sedov stage {Z is the particle charge). This finding may explain 
why the SNRs with the age more than a few thousand years are not prominent 
sources of very high energy 7-rays. The averaged spectrum of ultrarelativistic CR 
injected in the interstellar medium is close or somewhat steeper than E'"^. 

1. Introduction 

The dependent on energy diffusion coefficient D{E) determines the maxi- 
mum energy that particles can gain in the process of acceleration. The condition 
of efficient acceleration is D{E) < kushRsh, where i?sh is the radius and Ush is the 
velocity of spherical shock, = 0.1 in the free expansion stage and k = 0.04 in 
Sedov stage, e.g. [3]. The Bohm value Db = vr^/?) {v is the particle velocity, 
and Tg is the particle Larmor radius) that is a lower bound of the diffusion along 
the magnetic field gives -Emax = 1-7 x lO^'^Z {S^i/hq) eV at the end of the free 
expansion stage when particles reach the highest energy. Here the SN burst with 
kinetic energy of ejecta £ = ^silO^^ erg in the gas with density uq cm~^ and the 
interstellar magnetic field Bq = 5 x 10~^ G is considered. 

Analyzing the early stage of SNR evolution when the shock velocity is high, 
Ush ~ 10^ km s~^, it was found [2] that the CR streaming instability can be so 
strong that the amplified field 6B > 10~^ G far exceeds Bq. The maximum parti- 
cle energy increases accordingly. The CR streaming instability is less efficient as 
the shock velocity decreases with time and the nonlinear wave interactions reduce 
the level of turbulence at the late Sedov stage [7,10]. This leads to fast diffusion 
and decreases -Emax- The effect is aggravated by the possible wave damping on 
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the ion-neutral collisions [1,6]. In the present work, we consider the acceleration 
of CR and their streaming instability in a wide range of shock velocities. The 
level of magnetic field fluctuations is allowed to be arbitrarily large, and the rate 
of nonlinear wave interactions is assumed to correspond to the Kolmogorov non- 
hnearity. The coUisional dissipation is also taken into account. The task is to find 
the maximum energy of accelerated particles as a function of SNR age. 

2. Mciximum Energy of Accelerated Particles 

In the test particle approximation, the distribution of particles in momen- 
tum for high Mach number shocks has the canonical form f{p) ~ p^^. In the case 
of efficient acceleration, the action of CR pressure on the shock structure causes 
nonlinear modification of the shock that changes the shape of particle spectrum 
making it flatter at relativistic energies. So, we assume that the distribution at 
the shock is of the form fo{p) ~ where < a < 0.5, and value a — 0.3 

is used in the numerical estimates below. The normalization of function foip) is 
such that the integral N = An J dpp'^folp) gives the number density of CR. We 
assume that the CR pressure at the shock is some fraction ,^cr < 1 of the upstream 
momentum flux entering the shock front, so that = CcrP^gi^. The typical value 
of Ccr = 0.5 and the total compression ratio 7 were found for strongly modified 
shocks in [3]. 

The following steady-state equation determines the energy density W of 
the turbulence amplified by the streaming instability in the CR precursor up- 
stream of the shock: 

uvw = 2(rer - Ti - r,i)py. (1) 

Here the l.h.s. describes the advection of turbulence by supersonic gas flow. The 
terms on the r.h.s. of the equation describe respectively the wave amplification 
by CR, the linear damping of waves in background plasma and the nonlinear 
wave-wave interactions that limit the amplitude of turbulence. The Kolmogorov- 
type nonlinearity with the simplified expression Fni = {2CK)~^^'^Va_kA(> k) ^ 
0.05Vi^kA{> k) at Ck = 3.6 (as given by the numerical simulations [9]) is used in 
our calculations. Here A = 6B/Bq is the wave amplitude, k is the wave number. 
The wave-particle interaction is of resonant character and the resonance condition 
is kvg = 1. The equations for Fcr and for D used in our calculations generalize 
the standard equations [4] derived in the case of weak random field. The details 
of our consideration can be found in [8]. 

Fig.l illustrates the results of calculations of -Emax at the Sedov stage of 
SNR evolution at S = lO'""'^ erg in the warm interstellar gas with the temperature 
T = 8 X 10^ K, and the average density no = 0.4 cm . Three solid lines cor- 
respond to three cases of wave dissipation considered separately: nonlinear wave 
interactions; damping by ion-neutral collisions at constant gas density; damping 
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Fig. 1. The particle maximum momentum Pmax in units mc as a function of shock 
velocity Ugh- 

by ion-neutral collisions when the diffuse neutral gas restores its density after 
complete ionization by the radiation from SN burst. The maximum energy of 
protons accelerated by SN shocks at the early Sedov stage is close to 3 x 10^"^ eV 
(that exceeds the Bohm limit) and decreases to about 10^° eV at the end of the 
Sedov stage (that is less than the Bohm limit). In particular, the particle energy 
is less than 10^^ eV for t > 3 x 10^ yr and this may explain the absence of a 
TeV 7-ray signal from many SNRs [5] where 7-rays could be produced through 
7r° decays if sufficiently energetic CRs were present. 

The highest particle energy estimated as -Emax = 2 x 10^^Z(nsh/3 x 10^ km 
^~^yC,ciMll^n^/^ eV is reached at the end of the free expansion stage (Mgj is the 
mass of ejecta in solar masses). 

3. Average Spectrum of Injected Cosmic Rays 

Let us find the overall spectrum of CR injected to the interstellar space. 
It is not clear yet how the process of CR exit from the SNRs proceeds. Consider 
two cases: 

^4) All accelerated particles leaves the envelope without considerable adia- 
batic losses. The production of CR in the galactic disk is then Q = i^sn4vr / dtu^^^R'^^fo, 
where z/gn is the SN rate in the Galaxy. This gives the approximate scaling 
Q ~ J^snicr^P'^ (at p ^ mc) under the conditions that ^cr = const, and that 
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SNR expansion is adiabatic as in the case of Sedov solution. It is remarkable that 
the form of average source spectrum Q{p) is not sensitive to parameter a at a > 
(see also simulations [3]). 

B) The CR with maximum energies can not be confined near the shock and 
are leaving the precursor because of the Pmax(^) decrease. Other particles which 
are confined in the envelope experience adiabatic energy changes until the very end 
of Sedov stage and the exit from a SNR. The runaway particles has the average 
spectrum of the form whereas the trapped particles gives approximately 
p~^-^. The forms of these spectra are not sensitive to a at a > 0, and both 
populations have the comparable number densities. 

4. Conclusion 

The accounting for non-linear effects in the instability that accompanies 
the CR acceleration may simultaneously eliminate two difficulties of modern cos- 
mic ray astrophysics. It raises the maximum energy of accelerated particles in 
young SNR above the standard Bohm limit through the production of strong 
random magnetic fields and thus helps to explain the origin of galactic CR with 
energies up to VS^'^ Z eV. It also decreases the maximum energy of particles in the 
late Sedov stage of SNR evolution, which allows us to explain why these objects 
are not bright in very high energy 7-rays. It is remarkable that even with a fiat 
instantaneous particle spectrum at the shock typical for the strongly nonlinear 
regime of acceleration, the spectrum of ultrarelativistic particles injected into in- 
terstellar space and averaged over the age of SNR is close to or somewhat 
steeper. 

Acknowledgments. This work was supported by the RFBR grant. 

5. References 



1. Bell A.R. 1978, MNRAS 182, 147 

2. Bell A.R., Lucek S.G. 2001, MNRAS 321, 433 

3. Berezhko E.G., Yelshin V.K., Ksenofontov L.T. 1996, JETP 82, 1 

4. Berezinskii V.S. et al. 1990, Astrophysics of Cosmic Rays, (Amsterdam: North 
Holland) 

5. Buckley J.H. et al. 1998, A&A 329, 639 

6. Drury L. O'C, Duffy R, Kirk J.G. 1996, ApJ 309, 1002 

7. Fedorenko V.N. 1990, Preprint 1442, A.F.Ioffe Phys. Tech. Inst., Leningrad 

8. Ptuskin V.S., Zirakashvih V.N. 2003, |astro-p h/0302053| 

9. Verma M.K. et al. 1996, JGR 101, 21619 

10. Volk H.J., Zank L.A., Zank G.P. 1988, A&A 198, 274 



